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Abstract 

It is shown, that any sufficiently smooth periodic solution of the self-focusing Nonlin- 
ear Schrodinger equation can be approximated by finite-gap ones with an arbitrary 
small error. As a corollary an analogous result for the motion of closed curves in 
guided by the Filament equation is proved. This equation describes the dynamics 
of very thin filament vortices in a fluid. 



One of the basic questions of the finite-gap theory is the following: are the 
finite-gap solutions of a given equation sufficiently generic or they belong to 
some special subclass? To answer this question it is reasonable to check if 
arbitrary periodic in spatial variables solution can be approximated by finite- 
gap ones. 

The study of simplest examples shows, that the uniform approximation for 
all X and t is impossible, because to do it we have to keep all space and 
time frequencies simultaneously, and we have too many constraints on the 
deformations to fulfill all of them. 



Therefore it is natural to ask the following questions: 

1) Let us have a smooth periodic in spatial variables solution and an arbitrary 
compact domain U in the (x, t) space. Is it possible to approximate our solu- 
tion on U by finite-gap ones with arbitrary small error? The approximating 
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solutions are allowed to be non-periodic. We shall call such approximations 
local. 

2) Let us have a smooth periodic in spatial variables solution and an arbitrary 
compact domain U in the t space. Is it possible to approximate our solution 
by periodic in x finite-gap ones with the same x-periods for all x and t e U7 
We shall call such approximations periodic. 

Of course any periodic approximation is automatically local, but the charac- 
terization of periodic finite-gap solutions is usually sufficiently complicated, 
therefore the transition from local approximations to periodic ones may be 
rather nontrivial. 

In our text we construct periodic approximations for the self-focusing Non- 
linear Schrodinger equation (SfNLS) 

iqt + qxx + 2g^g = 0, where (1) 



q = q{x, t) is a complex- valued function of two real variables, and for the 
Filament equation 

^^^^k{x,t)b{s,t), where (2) 



7(s, t) is a t-dependent family of smooth curves in M^, s is the natural param- 

— * 

eter (i.e. \ds^{s,t)\ = 1), b{s,t), k{s,t) denote the binormal vector from the 
Prenet reper and the curvature function respectively. 

The first periodic approximation theorem was proved in 1975 by Marchenko 
and Ostrovskii [9] for the real Korteveg de Vries (KdV) equation. The method 
of [9] is based on the theory of conformal maps and it can be naturally extended 
to some soliton systems including the defocusing NLS. But for other systems 
the question is still open and the answer depends on the equation. For example, 
from results of Krichever [8] it is rather clear that any periodic Kadomthsev- 
Petviashvili II (KP II) solution allows finite-gap approximations, but it is 
likely that for KP I it is not so. 

Direct attempts to generalize the approach of [9] to SfNLS meet the following 
problems 

(1) The space of spectral curves corresponding to real periodic fif-gap KdV 
solutions is topologically M © (E+)^. But in the SfNLS case this space 
it the real part of some ramified covering of C^^^, and the structure of 
the Marchenko- Ostrovskii conformal map is essentially more complicated. 
To avoid a detailed study of the parameters space we use the method of 
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isopcriodic deformations suggested by M.Schmidt and author in [3]. 

(2) In the KdV case the characterization of admissible divisors is very simple: 
the n-th point of divisor is located at an arbitrary point of the n-th com- 
pact real oval. In the SfNLS case the characterization is less explicit (see 
below) and we have to describe how we vary the divisor after perturbing 
the spectral curve to preserve the admissibility. 

(3) Dubrovin equations for the real KdV are non-singular therefore a small 
change of parameters and starting point slightly affects the solution. But 
in the SfNLS case Dubrovin equations may have singularities, and the 
solutions may have branch points (nevertheless the corresponding SfNLS 
potential is smooth). Therefore we have to check that our variations of 
spectral data do not change the solution too much. To do it we introduce 
some new "symmetric" variables. 

Let us recall some basic facts from the SfNLS theory. The scattering trans- 
form for NLS was found in 1971 by Zakharov and Shabat [15]. Finite-gap NLS 
solutions were first constructed in 1976 by Its and Kotljarov [6]. The charac- 
terization of SfNLS admissible divisors as well as a proof that all solutions 
with reduction (6) are automatically nonsingular were obtained by Chered- 
nik [1]. Infinite-gap periodic problem for matrix operators including the NLS 
L-operator was studied by M.Schmidt [12]. A lot of useful information about 
the NLS theory including the Hamiltonian theory is contained in the book [2] 
by Faddeev and Takhtadjan. Finite-gap NLS theory is discussed in details in 
the article [11] by Previato. Effictivisation of low genus formulas by NLS was 
studied by Kamchatnov [7]. 

The zero- curvature representation for SfNLS reads as: 



a*(A..r.f) 



U{X,x,t)'i/{X,x,t), 



(9*(A.,r,t) 
di 



V{X,x,t)^{X,x,t), (3) 



dx 



where 'if{X,x,t) is a 2-component vector, 



^{X,x,t) 



iJi{X,x,t) 
ip2{X,x,t) 



(4) 



U{X,x,t), V{X,x,t) are the following 2x2 matrices: 



U{X,x,t) 



iX iq{x,t) 
ir{x,t) —iX 



,V{X,x,t) 



2XU{X,x,t) + 



iqr -q^ 



fx —iqr 



,(5) 



r(a;, t) = q{x, t). 



(6) 
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We shall assume that q{x, t) is periodic in x with the period 1 



q{x+l,t) = q{x,t). 



(7) 



We shall fix our attention on the spectral transform for a fixed t = to, therefore 
starting from this moment we shall omit t in our notations. 

The Bloch eigenfunction \E'(A,a;) is by definition the common eigenfunction of 
L = dx — U{X, x) and the shift operator 



Equation (8) defines ^(A, x) up to a constant factor. We fix it assuming 



The function p{X) is defined up to adding 27m, n G Z. It is called the quasi- 
momentum. 

To calculate the Bloch function we have to diagonalize the 2x2 monodromy 
matrix T(A), which is an entire function of A. The eigenfunctions of T(A) 
lie on a two-sheeted covering F of the A-plane. V is called spectral curve. 
Denote the permutation of sheets of F by a. detT(A) = 1, therefore p{'y) + 
p{a^) = 0(mod27r). (We shall denote points of F by 7 and the projection 
F — C by A = V'j). p{'-f) is a locally holomorphic multivalued function 
on F, dp = {dp{X) / d\)dX is a holomorphic differential on the finite part of 
F, a{dp) = —dp. p{'j) ~ ±A as A ^ 00 (as an asymptotic series), therefore 
F is compactified by 2 infinite points 00+, oo_, (700+ = oo_, (Too_ = 00+, 
p{'y) ~ ±A as 7 ^ oo± respectively. 

A point A e C is called reguleir if ^(7) 7^ 0(mod7r), where P7 = A and 
irregular otherwise. Let Afc be an irregular point. The Tailor expansion of 
trT(A) reads as: trr(A) = ±2 + T^{X - A^)" + . . .. Let us call n the order 
of the point Afc, n = ordg(Afc). {q means that the order is defined in terms of 
the quasimomentum function). A^ is a branch point of F if n is odd and a 
double point of F if n is even. A branch point is called simple if n = 1. If 
the opposite is not stated exphcitly we have one point of F over each double 
point. To have an uniform representation for our equations wc shall treat an 
irregular point of order n as the result of fusion n simple branch point. 

If ^(A,a;) is a Bloch solution of (8), then 



L*(A, x) = 0, *(A, X + 1) = e^f(^)*(A, x). 



(8) 



$(A, 0) = 1, where ^{\,x) = ipi{\,x) + ip2{X,x). 



(9) 



^+(A,a;) 



V'2(A,x) 
-V5i(A,x) 



(10) 
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is also a Bloch solution of (8) with the quasimomentump"'"(A) = —p{X). There- 
fore r has the following antiholomorphic involutions 7 — (77 and 7 — 7. (We 
assume that ^(7) = ^(7), p(cr7) = — ^(7)). For real A T(A) is an unitary ma- 
trix, p(A) G K, and V has no real branch points (but may have real double 
points) . 

Let {Ek}, {E^} be the lists of all irregular points, the index k takes all integer 
values. We assume that 

(1) ImEk > 0. 

(2) E+ = e;. 

(3) If ordq(Afc) = n the point A^ has exactly n entries in our lists. For example 
if ordg(Afc) — 4 and A/; G K then we have exactly 2 integers /ci, /c2 such 
th^t Ek,^E,,^E+^E+. 

Lemma 1 It is possible to enumerate the irregular points so, that for suffi- 
ciently large \k\ 

(1) Ek^ (TTSgn/c) • - I-^{q) + o (i) where Ii{q) = Jq q{x)q{x)dx. 

(2) Ek — E^ — > faster than any degree of k'^ as \k\ — > 00 (We assume q{x) 
to he smooth). 

The next important object for us is the set of zeroes of the quasimomen- 
tum differential dp. They are invariant under the involution a therefore we 
shall consider their projections to the A-plane instead. Denote them by 
where k takes all integer values. As above we use the following agreement 

(1) If Afc is a regular points it has n entries in the list {a^} where n is the 
order of zero of dp at one sheet. 

(2) If Afc is an irregular points it has 
[ ] denotes the integer part. 

Lemma 2 It is possible to enumerate the points so, that for sufficiently 
large \k\ 

(1) lmafc = 0. 

(2) ak = ReEk + o{lmEk). 

Let us define now the "second part" of the spectral data - the divisor of 
poles of the Bloch function. 

Let ^(7,0;) denote a Bloch eigenfunction of L with some non-singular lo- 
cally holomorphic normalisation (of course ^(7,0;) = ^(7, a;)/$(7, 0) where 
$(7,0;) = V^i(7,x) + 1^2(1, x). Consider the Wronskian of the Bloch functions 
1^(7) = ilji{'j,x)ip2{c!"l,x) — ip2{l,x)ipi{aj,x). It is defined up to a non-zero 
holomorphic multiplier and does not wanish at regular points. Let A^ be an 



entries to the list {ak} where 
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irregular point. We have W{\) = ±w™(A — Afc)™'/^(1 + o(l)) where m is even 
if Afc is a double point and odd if is a branch point, m > 0. Denote m 
by ordb(Afc). It is easy to check that ord{,(Afc) < ordg(Afc). A double point A*; 
is called removable if ord6(Afc) = 0. It is well-known that removable double 
points can be treated as regular points and we can forget about them. 

The divisor of Bloch function zeroes is a list of points of F {'jk{x)} 
where k takes all integer values such that each zero of $(7,0;) generates / 
entries to this list if I is the multiplicity of it and each irregular point generates 
(ord^ — ord6)/2 entries. The divisor of Bloch function poles {7^} coinsides 
with the divisor of Bloch function zeroes taken at the point x — 0. 

Lemma 3 (1) The spectral curve F has only finite number of non-removable 

double points and degenerate branch points. 
(2) All real double points are removable. 

Lemma 4 It is possible to enumerate the points 7^ so, that for sufficiently 
large \k\ V-fk ^ReEk + 0{lmEk). 

It is well-known, that the spectral curve and the divisor of poles completely 
define the potential q{x). To reconstruct the potential we can use Dubrovin 
equations 



fe=— 00 



u,{x), (11) 



where 



^ji^) = yM^) - Ej){\j{x) - E+) n x.(^\_x '(12) 

k^j i\ ) k\ ) 

and the reconstruction formula 

q{x)= E P^-A.(a;) + E ^.(^)- (13) 

fc=— 00 \ / j=—oo 

The infinite sums and products in the formulas above perfectly converge. 

Let us recall the characterization of divisors corresponding to operators with 
the reduction (6). Consider the following 1-form on F: Q{'y,x) — uj{'y,x)dX, 
where 

/ ^ +ilj2{-f,x)){'ilj2{a-f,x) -iJi{a^/,x)) 

UJ['-f,x) — = = = = . (14j 

A{l,x)iJ2{(Tl,x) -^2(7,a;)^i(cr7,a;) 
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Denote by U{R) be the domain |A| < in F. Consider the following function 
u}R{'y,x) — ou{'y,x) \/(A — E'^) (A — E'j^). Denote by ^(a;^, x) the divisor 

fc:|Efe|<i? 

of zeroes of u!ii{'y, x) and by D{u!, x) the limit of D{u!r, x) as — > oo. 

Lemma 5 D{lj,x) coincide with the set {'Jk{^),Jk{-'^)}, where {■Jki^)} is the 
divisor of Bloch function zeroes. 

Let 5j be an arbitratry collection of pairwise distinct real points such that for 
sufficiently large \j\ Sj — ReEj. 

Lemma 6 (1) The form Q reads as Q = [1 — x)] dX 

~ ^( A - E,) (X-El) t}j ^{X - Ek) (A - O ' 

where Kj{x) are some real functions of x. 
(2) |k(7,x)| < 1 for all X eR, J e M. 



functions 



In particular K,[5j,x)\ < 1 for all j. It gives us the following estimate on the 



\Kj{x)\ < 



k^j ^3 - 



<\5,-E,\Cr,{16) 



where Cr is a positive constant, depending only on the spectral curve. In 
particular, if we have a removable double point Ek — E^ — 5k, then Kk{x) = 0. 

Lemma 7 Let ^^(0) be a collection of real numbers such, that \k{j,0)\ < 1 
for all 7 G r, where fi;(7, 0) is defined by (15), D{uj,0) be the corresponding 
divisor, {7j(0)} be any set of points such that D{lj,0) — {7fe(0), 7^(0)}. Then 
the corresponding operator L-operator satisfy (6), and the potential q{x) is 
nonsingular. 

For us the following definition will be convenient: potential q{x) is called 
finite-gap if Ej — Ej' for all \j\ > Jq. Then all points Ej — Ej' are re- 
movable double points, aj = 7j(a;) = Ej for all \j\ > Jq and F has only 
finite number of branch points and non-removable double points. Finite-gap 
solutions of soliton equations were first introduced by Novikov in 1974 for 
KdV [10]. The corresponding solutions can be written explicitly in terms of 
Riemann ^-functions. 

The first step of the approximation procedure is to construct a finite-gap 
deformation of F generating solutions with the same period. To do it we need 
the following lemma proved by M.Schmidt and the author in [3]. 
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Lemma 8 Let ak ^ M. be the projection of a zero of the quasimomentum such, 
that «fc 7^ aj for j ^ k, ^ Ej, ^ Ej~, for all j. Consider the following 
system of ODE's, associated with the point a^; 



dE^ 
dr 



^ ^ 

Ea -ak dr 



Ck{r) 



Ej - ak ' 



daj 
97 



Cfc(T) 



Oik 



for j ^ k, 



(17) 



dak 
dr 



Ck{r) 




where Ckir) is an arbitrary real function ofr. 

Denote by T{r) the solution of (17) with the initial value r(0) = T, where the 
spectral curve V corresponds to a periodic with the period 1 potential q{x) (of 
course this solution is defined only in some neighborhood of zero U{0)). Then 
for all T G f/(0) the curve T{r) generates periodic with the period 1 potentials 
(the x-quasifrequencies of the potentials do not depend on the divisor). 

Let I A; I be sufficiently large. Then using this deformation we can merge the 
pair Ek, E^ to a removable double point, and the corresponding shift of all 
points Ej, Ej~, Uj with j ^ k is of order o{lmEk). Therefore applying this 
deformation to all k such that \k\ > K, where K > is a sufficiently large 
integer, we obtain a finite-gap spectral curve Tk (it is almost evident that the 
superposition of infinitely many deformations perfectly converges). 

Lemma 9 For any e > there exists a K such, that 

(1) \Ej — Ej\ < e, \Ej' — Ej~\ < e, \aj — aj\ < e, for all j. 

(2) I lra{Ej - Ej)\ < e\ Im Ej\, for all j such, that \ j\ < K. 

where Ej, E^ , aj are the branch points of the curve and the quasimomen- 
tum zeroes respectively. 

We have constructed a family of finite-gap curves approximating the curve 
r. Let us discuss now the admissible divisors. 

Lemma 10 There exists a pair of positive integer constants Ki, K2 such, 
that for all K > K2 the points of any admissible divisor 7^ on can be 
enumerated so, that 

(1) For all k such that \k\ < Ki \\k\ <Ki + 1/10. 

(2) For allk such that \k\ > Ki \{Xk-Sk)\ < Im-Bfc, and |4| > -f^i + l-l/lO. 



8 



The proof follows from the characterization of admissible divisors given by 
Lemmas 5-7. 

Equations (11)- (12) have singularities at the right-hand side. To simphfy the 
structure of Dubrovin equation it is convenient to introduce the following new 
variables: 

(1) Sfe(x), qk{x), 1 < k < 2Ki + 1 - the first 2Ki + 1 expansion coefficients 
at oo of the function 

^ J(X- Ek)(X- E+) 
Eh.x) = (l + .(,,.))nVi_ji_^ (18) 

^(A, x) = ±\ l + }_^ -—- + 2^ -—- as 7 ^ ±oo. (19) 

\ k>0 J k>0 

(2) \k{x) = Xk{x)-Sk, \k\ > Ki. 

(3) M^) = ^{Xk{x)-Ek){Xk{x)-E^), \k\ > K^. 

These variables are dependent. Denote this set of variables by S. 
Consider the following norm: 

ll'^lln = .1 E {\Sk?+W)+ E (20) 
V \k\<Ki \k\>Ki 

This norm is bounded on the space of admissible divisors for any positive n. 

Lemma 11 For any sufficiently large n there exists a constant CniT) such, 
that for any admissible pair Si {x) , S2 {x) of solutions of Dubrovin equations 
we have the following estimate 

<Cn{V)\\Si{x)-S2{x)\\n. (21) 

n 

It is easy to check, that for any ei > there exists a constant Ky,{n) such that 
||5||(f)<6i, where = ./ E \k\- {\\k? + W) ■ (22) 

V \k\>K3(n) 

We need also the following semi-norm 

ll'5||« = J E (\sk? + \qk?)+ E \k\-{\Xk? + H?). (23) 

V |fc|<A'i Ki<\k\<K3{n) 



-iSi{x) - S,{x)) 
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It is evident, that + \\S\\!^^ and < \\S 



n- 



Combining (21) and (22) we obtain the following estimate 



S,{x) - S,{x)\U < e^"(^)N (||5,(0) - 52(0)||« + e,) . 



(24) 



Therefore in approximate calculations we can truncate the Dubrovin system 
to a finite- dimensional one, removing the variables Xk{x), i'kix), with the 
\k\ > Ki{n). For the truncated system small variations of the curve and of 
the starting point result in small variations of the solution. To complete the 
proof it is sufficient to check, that choosing K sufficiently large we can make 
the admissible variation of divisor arbitrary small. But it follows from the 
characterization of admissible divisors presented above. 

These arguments can be applied also for the Dubrovin equations, describing 
the ti -evolution of the divisor, where ti denotes the l-s time from the NLS 
hierarchy (in these notations t = ti). Taking into account that the first k x- 
derivatives of q{x) are continuous functionals in the norm || ||n for sufficiently 
large n we obtain: 

Theorem 1 Let q{x, t) be an arbitrary SfNLS solution with smooth x-periodic 
Cauchy data q{x,0) = qo{x). Then for any e > 0, N > and T > there 
exists a finite- gap SfNLS solution q^{x,t) such, that 



At the end let us say a few words how to prove an analogous theorem for the 
Filament equation. Equivalence between the SfNLS and the Filament equation 
is given by the Hasimoto map [5] . (^-functional solutions of Filament equations 
were studied by Sym in [13]). 



where k{s,t), K.{s,t) are the curvature and the torsion functions respectively. 
In [4] it was shown, that equations (17) (except one corresponding to ao = 0) 
preserve the periodicity in s of the Filament equation solution. Therefore the 
technique developed above can be applied without changes. 

Acknowledgments. The author is grateful to Martin Schmidt for explaining 
some results from the periodic NLS theory. 




for all X eR, \t\<T, n< N. (25) 




(26) 
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